Current-voltage characteristics of a graphene nanoribbon field-effect transistor 
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We present an analytical device model for a field-effect transistor based on a heterostructure which 
consists of an array of nanoribbons clad between the highly conducting substrate (the back gate) 
and the top gate controlling the source-drain current. The equations of the model of a graphene 
nanoribbon field-effect transistor (GNR-FET) include the Poisson equation in the weak nonlocality 
approximation. Using this model, we find explicit analytical formulas for the spatial distributions of 
the electric potential along the channel and for the GNR-FET current-voltage characteristics (the 
dependences of the source-drain current on the drain voltages as well as on the back gate and top 
gate voltages) for different geometric parameters of the device. It is shown that the shortening of the 
top gate can result in a substantial modification of the GNR-FET current-voltage characteristics. 
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I. INTRODUCTION 

Graphene, i.e., a monolayer of carbon atoms forming 
a dense honeycomb two-dimensional (2D) crystal struc- 
ture is considered as a promisingcandidate for future 
micro- and nanoelectronics [J, [3,11, 0, S] • The features 
of the electron and hole energy spectra in graphene pro- 
vide the exceptional properties of graphene-based het- 
erostructures and devices, in particular field-effect tran- 
sistors [3, Fabrication and experimental studies as 
well as a theoretical model of graphene field-effect tran- 
sistor (G-FET) was reported recently in, for instance, 
Ref. [□]. The operation of G-FETs is accompanied 
by the formation of the lateral n-p-n (or p-n-p) junc- 
tion under the controlling (top) gate and the pertinent 
energy barrier. The current through this barrier can 
be associated by both thermionic and tunneling pro- 
cesses H, 0, B S E3 ■ The utilization of the patterned 
graphene which constitutes an array of sufficiently nar- 
row graphene strips (nanoribbons) provides an opportu- 
nity to engineer the band structure to achieve the optimal 
device parameters. In particular, properly choosing the 
width of the nanoribbons, one can fabricate the graphene 
structures with relatively wide band gap but rather high 
electron (hole) mobility 0, [H|. 

In this paper, we present an analytical device model 
for a graphene-nanoribbon FET (GNR-FET) and obtain 
the device characteristics. The GNR-FET under con- 
sideration is based on a patterned graphene layer which 
constitutes a dense array of parallel nanoribbons of width 
d with the spacing between the nanoribbons d s <§; d. The 
nanoribbon edges are connected to the conducting pads 
serving as the transistor source and drain. A highly con- 
ducting substrate plays the role of the back gate, whereas 
the top gate serves to control the source-drain current. 
The device structure is schematically shown in Fig. 1. 
Using the developed model, we calculate the potential 
distributions in the GNR-FET as a function of the back 
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FIG. 1: Schematic side (a) and top (b) views of a GNR-FET 
structure. 



gate, top gate, and drain voltages, Vb, V g , and Vd (reck- 
oned from the potential of the source contact), respec- 
tively, and the GNR-FET dc characteristics. This cor- 
responds to a GNR-FET in the common-source circuit. 
The case of common drain will briefly be discussed as 
well. For the sake of definiteness, the back gate voltage 
Vb and the top gate voltage V g are assumed to be, re- 
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spectively, positive and negative (Vb > and V g < 0) 
with respect to the potential of the source contact, so we 
consider GNR-FETs with the channel of n-type. 

The paper is organized as follows. In Sec. II, the GNR- 
FET model is considered and the main equation gov- 
erning the spatial distributions of the electric potential 
in the active region of the transistor channel is given. 
This equation is a consequence of the two-dimensional 
Poisson equation in the weak nonlocality approxima- 
tion [l(| EH- Section III deals with the derivation 
of the relations between the electron and hole densities 
in the channel and the electric potential and the con- 
sequent calculations of the potential spatial distributions 
depending on the drain and gate (back gate and top gate) 
voltages. In Sec. IV, using the obtained formulas for the 
electric potential and the height of the barrier for elec- 
trons at the minimum of the potential, we obtain analyt- 
ical expressions for the source-drain current as a function 
of the drain and gate voltages. In Sec. V, we discuss the 
specific of GNR-FET operation in the circuits with com- 
mon drain, the possible role of holes, and the limitations 
of the model used. In Sec. VI, we draw the main results. 



band nanoribbons. Notwithstanding this, it is suggested 
that the channel sections adjacent to the source and drain 
contacts are highly conducting, so that these section are 
equipotential. The potentials are equal to the potentials 
of the source and drain contacts, i.e., equal to ip — and 
ip = Vd, respectively. 

We shall consider the GNR-FET region under the 
top gate (the device active region) defined as follows: 
—L g < x < L g /2, —Wb < z < W g , where L g is the 
length of the top gate, and Wb and W g are the thick- 
nesses of the layers between graphcne and the back and 
top gates, respectively. Here, the axis x is directed along 
the nanoribbons, whereas the axis z is directed perpen- 
dicular to the nanoribbons and gate planes. 

To find the potential distribution along the channel, 
we use the following equation: 

{W h + W g )#<p V -V b if-Vg 

3 dx 2 W b W g 33 1 +> 

(2) 

with the boundary conditions 



II. GNR-FET MODEL 

Graphene strips (nanoribbons) exhibit the energy spec- 
trum with a gap between the valence and conduction 
bands depending on the nanoribbon width d: 

£j,n = ±« vV + (Trh/d) 2 n 2 . (1) 

Here v ~ 10 8 cm/s is the characteristic velocity of the 
electron (upper sign) and hole (lower sign) spectra, p is 
the momentum in along the nanoribbon, H is the reduced 
Planck constant, and n = 1,2, 3, ... is the subband index. 
The quantization corresponding to Eq. (1) of the elec- 
tron and hole energy spectra in nanoribbons due to the 
electron and hole confinement in one of the lateral direc- 
tions results in the appearance of the band gap between 
the valence and conduction bands and in a specific the 
density of states (DOS) as a function of the energy. 

The electron density in different parts of the channel 
and, therefore, the degree of the electron gas degener- 
ation essentially depend on the back gate voltage. In 
contrast to graphene with zero energy gap, the electron 
(hole) gas in relatively narrow nanoribbons with quan- 
tized energy spectrum and the pertinent energy gap be- 
tween the subbands in the valence and conduction band, 
can become degenerate at fairly high back gate voltages. 
In the following, we restrict ourselves by the considera- 
tion of the GNR-FET operation under the condition that 
the electron gas in the channel is nondegenerate. Such a 
consideration is valid in the voltage range which is suffi- 
ciently wide in terms of practical applications. Thus, the 
back gate voltage is assumed to be not excessively high, 
so that the electron density is moderate, the electron gas 
in the channel is nondegenerate, and the electrons oc- 
cupy only the lowest [n = 1) subband in the conduction 



<f\x=-L g /2 = 0, <P\x=L a /2 = Vd- (3) 

Here, £_ and S+ are the electron and hole sheet den- 
sities in the channel, e is the electron charge, and ae 
is the dielectric constant of the material separating the 
channel from the gates. Equation (2), which governs the 
electric potential, tp = <p(x) = ip(x,0), in the channel 
is a consequence of the two-dimensional Poisson equa- 
tion for the electric potential ip = ip(x, z) for the active 
region under consideration in the weak nonlocality ap- 
proximation [l2l . Il3j |. This equation provides solutions, 
which can be obtained from the exact solution of the 
two-dimensional Poisson by the expansion in powers of 
the parameter <5 = [(W b 3 + W^)/45(W b + W g )C 2 }, which 
is much smaller than unity in the situation under consid- 
eration. Here J is the characteristic scale of the lateral 
inhomogeneities (see Sec. V). The lowest approximation 
in such an expansion corresponds to the so-called grad- 
ual channel approximation proposed by W. Shockley, in 
which the first term in the left-hand side of Eq. (2) is 
neglected [3, [T^], so that the relationship between the 
potential in the channel and the electron and hole charge 
becomes local. Thus, Eq. (2) can be used when <5 < 1 and, 
hence, when C > max{Wb, Wg}. The harnessing of the 
approximation under consideration makes, in particular, 
possible to study essentially nonuniform potential dis- 
tributions in the GNR-FET channel and the short-gate 
effects analytically. The weak nonlocality approximation 
was effectively used previously in some previous papers 
(see, for instance, Refs. [ICt Il2i. Il3l. Il6||). 

Since the spacing, d s between the nanoribbons is small, 
we disregard a small scale nonuniformity of the electric 
potential distribution in the in-plane direction y. 



3 



III. POTENTIAL DISTRIBUTIONS 

The application of negative top gate voltage leads to 
the formation of a potential barrier in the channel under 
this gate. This barrier determined by the gate voltage 
controls the source-drain current and, hence, is responsi- 
ble of the device operation as a transistor. 

A. Electron and hole densities 

Considering that the electron and hole gases are non- 
degenerate, and, hence, the electron and hole distribution 
functions in the subbands with n = 1 near the source and 
drain is given by 



/ ±e<p ±e F - y/v 2 p 2 



- A 2 /4 



where s F is the Fermi energy reckoned from the middle of 
the energy band gap and A = Eqj — e^i = 27r vh/d is the 
value of the energy band gap, the electric potential, the 
electron densities and the Fermi energies in the source 
and drain regions (marked by superscripts "s" and "<i" , 
respectively) are related to each other as 



eV b 



16 V^k^T\W b 



h v 



(7) 



sinh ( ' 



A~2eV d \ 



33 e{V b -V d ) / d \/Hv 
160F VAkFT \WbJ\e 2 

When Vb — V g — Vd — 0, one obtains e s F 



-d 
- F 



(8) 



0, 



so that the electron and hole densities are equal to their 
thermal value ujh 



2y/Ak B T 
y^nhdi 



exp 



2k B T 



(9) 



For the effective operation of a GNR-FET with the n- 
type channel, the electron densities in the sections of the 
channel should be reasonably high (to provide necessary 
conductivity of these sections). This requires the appli- 
cation of a sufficiently high back gate voltage. If 



exp 



A 



2k B T 
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Here, £ m = A/2k B T and Ki(£) is the modified Bessel 
function. Taking into account the asymptotic behav- 
ior of the Bessel function at £ ^> 1, namely Ki(^) ~ 
(7r/2£) 1/2 exp(-£), for A > 2k B T, we obtain 



^ y/nhdv 
2y/Ak B T 



exp 



etp 



k B T 



i/irhdv 



■ exp 



±- 



e(ip 



2k B Tj : 
V D ) 



k B T 



2k B T 



(6) 

Considering that when ip = and ip = V d , £1 — ^+ o = 
tsiVbl^eWb and E£ )0 - Ef )0 = ae (V b - V d ) /4tt eW b , re- 
spectively, where the quantities with the index "0" are 
the electron and hole densities in the immediate vicinity 
of the source and drain contacts, we arrive at the follow- 
ing equation: 



sinh 



k B T) 



exp 



A 



2k B T 



where 



V F = 



VAfc^T 



8^ rw b \f e 2 

se \ d J \hv 



(11) 



(4) as follows from Eqs. (7) and (8), s s F ,e F > k B T, and 



2 B \V F 

* + k B TlJY^t 
2 B \ V F 



(12) 



At Vb and (V b — V d ) satisfying inequality (10), the electron 
gas is nondegenerate and the electron density markedly 
exceeds its thermal value as well as the density of holes. 
In the case under consideration, the electron density in 
the pertinent sections of the channel, as follows from 
Eqs. (6) and (12), are given by 



33 V b 

exp 

AireWb 

33 (Vb - Vd) 
AireWb 



eip 
k^T 



exp 



e(<p - V d ) 



k B T 



(13) 



Since at Vb > 0, the electron density markedly exceeds 
the hole density, we can neglect E^ in the right-hand side 
of Eq. (2). The case when the hole gas under the top 
gate becomes essential (at very high top gate voltages) 
will be discussed in the following (in Section V) . Consid- 
ering the relationships between the electron density and 
the electric potential, from Eq. (2) one can arrive at the 
following equations governing the potential distribution 
in the active region: 
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FIG. 4: Spatial distributions of the potential in GNR-FETs 
with different gate lengths. 
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FIG. 3: Spatial distributions of the potential at different drain 
voltages. 



B. Potential distributions at low top-gate voltages 

When the top gate is negative (V g < 0) and its abso- 
lute value \V g \ is sufficiently small, the modulus of the 
potential \tp\ can be not that large. In this case, one can 
linearize Eqs. (14) and (15) and present these equations 
in the following form: 



d 2 V 
dx 2 



W b W q 



eV b 



(W b + Wg) k B T 



3V„ 



{W b + W B )W g ' 



(16) 



ci 2 cp 

Jx 2 ~ W b W,_ 



3 = _ 3(V b /W b + V g /W g ) 



(W b + Wg) 



(W b + Wg)W b 



V b ( eip 
exp 



d 2 y 3 3(V b /W b + Vg/Wg) 

dX 2 WbWg^ (Wb + W g ) 



dx 2 



WbWn 



w n 



e(V b - V d ) 



(W b + Wg) k B T 
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HVg/w g ~ v d /Wb) 

(Wb + Wg) 



(17) 



(14) Equations (16) and (17) should be solved taking into 

account Eq. (3) and the matching of (p and dtp/dx at a 



certain point x 



(such that —Lg/2 < x m < L g /2), 



at which tp reaches a minimum. Such a minimum defi- 
nitely exist when V g < 0. At V d = 0, Eqs. (16) and (17) 
yield 



(V b - V d ) 



(W b + W g 



W b 



■ exp 



e(<p - V d ) 



knT 



(15) 



The exponential dependences in the right-hand sides of 
Eqs. (14) and (15) are valid provided that inequality (10) 
is satisfied, in particular, if the electron gas is nondegen- 
erate. The threshold value of the back gate voltage, at 
which the degeneration of the electron gas occurs, can be 
estimated as Vb — Vf- 



VgW b 
(W b + Wg) 



where 



A = 



cosh(a;/A) 
cosh(L ff /2A) 



W n 



eV b 



(W b + W g ) k B T 



(18) 
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is the effective screening length and A = y/WbW g /3. At 
Vd = 0, as follows from Eq. (18), the function ip exhibits 
a minimum p — ip m0 at x — with 



<PmO 



VgW b 

(W b + Wg 



1 - 



1 



cosh(L g /2A) 



1 + 



W a 



eV b 



(W b + W g ) k B T 



Vg W b k B T 

V b Wg e ' 



(19) 



Here, we have taken into account that normally Vb S> 
k B T/e and L g > A (with A < A). 

At reasonable values of the drain voltages Vi, the vari- 
ation of p m — ip m0 is equal to Vd/2 cosh(L 9 /2A) <C Vd, 
i.e., is very small due to cosh(L s /2A) ;» 1 and can be 
disregarded, so that 



Pm — <PmO — 



V g W b k B T 
V b W q e ■ 



(20) 



Comparing e\ip m \ given by Eq. (20) with k B T, we find 
that Eqs. (16) - (19) are valid when \V g \ < V b W g /W b 



C. Potential distributions at high top-gate voltages 

At high back gate voltages, the quantity A playing the 
role of the screening length is rather small. In this case, 
the length of the regions near the points x = ±L g /2, in 
which the potential changes from <p = to \<p\ > k B T 
and from (p = Vd to \ip — Vd\ > k B T, is small in compar- 
ison with the top gate length L g . In such short regions 
near x — ±L g /2, the potential distribution can still be 
describe by Eq. (18). However, in a significant part of 
the active region the electron charge, which is in such 
a situation exponentially small, can be disregarded and 
the last (exponential) terms in Eqs. (14) and (15) can 
be omitted. Taking this into account, at high top gate 
voltages, Eqs. (14) and (15) in the central region can be 
presented in the following form: 



dx 2 



W b W c 



■p = — 



H Vg/Wg + Vb/W b ) 

(W b + Wg) 



(21) 



Solving Eqs. (21) still using boundary conditions given 
by Eq. (3), we arrive at 



Vg + Vb 



Wg\ Wb_ 



Wb) (Wb + Wg) 



1 - 



cosh(a;/A) 



cosh(L s /2A) 



+ V d 



sinh[(a; + i 9 /2)/A] 



sinh(L ff /A) 

At Vd = 0, ip exhibits a minimum at x — and 



(22) 



PmO ^ [ Vg + Vb 



Wg\ W^ 



'Wb) (Wb + Wg 



cosh(L 9 /2A) 



(23) 

In the limit L g ^> A, the potential in the center of the 
channel given by Eq. (22) coincides with that obtained 
in the gradual channel approximation. However, in real 
GNR-FET structures, L g can be comparable with A. 
In this case, the terms dependent on L g can be impor- 
tant (in contrast to the case of low gate voltages consid- 
ered in the previous subsection). At reasonable values of 
the drain voltage Vd (sufficiently small compared to Vb), 
Eq. (22) yields 



V >V W A Wb 

V9+Vb W b ) (Wb + Wg) 



v d 



1 - 



1 



cosh(L 9 /2A) 



2cosh(L s /2A)' 



(24) 



Figures 2 and 3 show examples of the spatial distri- 
butions (along the channel, i.e., in the x-direction) of 
the electric potential in the active region (under the top 
gate) calculated for a GNR-FET with W b = 100 nm, 
W g = 30 nm, and L g = 300 nm at the back gate volt- 
age Vb = 2.0 V, assuming different values of the top gate 
voltage V g and the drain voltage Vd- As seen from Fig. 2, 
the potential distribution markedly sags and the height 
of the barrier for electrons near the center increases with 
increasing absolute value of the top gate potential \V g \. 
Figure 3 demonstrates, in particular, that in the GNR- 
FET with chosen parameters the minimum value of the 
potential ip m and, hence, the height of the barrier — ep m 
for the electrons propagating from the source are vir- 
tually insensitive to the drain voltage. In contrast, the 
potential barrier e(Vd — p m ) for the electrons propagat- 
ing from the drain increases linearly with increasing Vd- 
Figure 4 shows the spatial distributions of the electric po- 
tential calculated for the GNR-FETs with different gate 
lengths L g . 

The approach used in this subsection is valid when 
e\p m \ > k B T, i.e., when \V g \ - V b W g /W b > k B T(W b + 
Wg)/eW b ~ k B T/e. 

In the limit L g 3> A > A the equations for p m ob- 
tained above (in particular Eq. (24)) coincide with those 
obtained using the gradual channel approximation. 



IV. CALCULATION AND ANALYSIS OF 
CURRENT- VOLTAGE CHARACTERISTICS 

Considering that the source-drain current is deter- 
mined by the electrons overcoming the potential barrier 
under the top gate, one can use the following formula for 
the density of this current (per unit length): 



J = 
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FIG. 5: The source-drain current density versus drain voltage 
dependencies at fixed back gate voltage (Vb = 2.0 V) and 
different top gate voltages V g . 



Here. 



ut p.O 2 



y/v 2 P 2 + A 2 /A 



is the velocity of the electron with momentum p in the 
lowest subband of the nanoribbon conduction band and 
p s m and p^ are the momenta of the electrons with the 
energies e\ip m \ and e(|<p m | + Vd), respectively. We have 
taken into account that the nanoribbon array is dense: 
d s <C d. Integrating in Eq. (26), we arrive at 



2tt 3 / 2 W b 



k B T 



■ exp 



etpr, 



k B T 



V b -{V b -V d )exp - 



eV d 
k B T 



(26) 



with tp m given by Eq. (20) at low top gate voltages and 
by Eq. (24) at moderate and high voltages, respectively. 
Equation (26) corresponds to the average thermal elec- 
tron velocity vt = v^J Ak B T / tt A. One can see that the 
thermal electron velocity as well as the pre-exponcntial 
factor in the formula for the current depend on the band 
gap. This is due to the dependence on the gap of the 
electron dispersion and the electron velocity. 



A. Current versus drain voltage 

The dependence of the source-drain current on the 
drain voltage is associated with the dependence of tp m 
on this voltage given by Eqs. (20) and (24) and the volt- 
age dependence of the last factor in the right-hand side 
of Eq. (26). At low top gate voltages (\V g \ < V b W g /W b 
), using Eqs. (20) and (26), we obtain 



J oc 



(27) 



FIG. 6: The same as in Fig. 5 but for V b = 3.0 V. 



i.e. 



if V d < k B T/e , and 



Jcxl/ rf 



J = const 



(28) 



(29) 



if k B T/e « V d < V b W g /W b . 

At moderate and high top gate voltages when (p m is 
given by Eq. (24), we arrive at the following dependence 



eV d 



J tx exp 



V b -{V b -V d )exp[- 



2k B Tcosh(L g /2A) 
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(30) 



At V d < k B T/e, Eq. (30) yields the same linear depen- 
dence on the drain voltage as that described by Eq. (28). 
When V d 3> k B T/e, one obtains 



J oc exp 



eV d 



2fc B Tcosh(£ g /2A) 



(31) 



The latter dependence resembles the dependence given 
by Eq. (29), but with the substitution of A by A. Since 
A <C A, the source-drain current at moderate and high 
top gate voltages is a less steeper function of the drain 
voltage than that at low top gate voltages (compare 
Eqs. (29) and (31)). 

The J versus V d dependences for a GNR-FET with 
L g = 300 nm at the back gate voltages V g = 2 V and 
V g = 3 V calculated for the different top gate voltages V g 
are shown in Figs 5 and 6. Here, as in the previous and 
consequent figures, we assumed that A = 0.4 eV, ae = 4, 
W b = 100 nm, W g = 30 nm, and T = 300 K. Since the 
gate voltages were set to be relatively high (V b , \ V g \ ^> 
k B T/e ~ 0.025V), Eqs. (24) and (26) were used for the 
calculations. 

As seen from Figs. 5 and 6, the source-drain current 
as a function of the drain voltage in a GNR-FET with 
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FIG. 7: The source-drain current density as a function of the 
top gate voltage at different back gate voltages. 



relatively long gate (L g — 300 nm) exhibits saturation 
starting rather low drain voltages. This is a consequence 
of very weak sensitivity of the potential barrier for the 
electrons propagating from the source to the drain volt- 
age. 



B. Current versus gate voltages 

As follows from Eqs. (20), (24), and (26), the source- 
drain current as a function of the top gate voltage is 
described by the following relations: 



(32) 



at \V g \<V b W g /W b , and 



J oc exp 



eV a 



k B T (W b + W g 



1 - 



cosh(L 9 /2A) 



(33) 



at \V g \ > V b W g /W b . 

Figures 5 and 6 show the transformation of the J vs Vd 
dependences with varying top gate voltage in the range 
of the latter \V g \ > V b W g /W b . These figures confirm a 
strong sensitivity of the source-drain current to the gate 
voltages. As follows from Eqs. (32) and (33), the depen- 
dence of the source-drain current on the top gate volt- 
age is much steeper in the range high top gate voltages 
\Vg\ > VbWg/Wb than at \V g \ < V b W g /W b , i.e., when 
the central region of the channel becomes essentially de- 
pleted. Figure 7 shows the J versus V g characteristics 
obtained using Eq. (26) in which tp m was determined by 
Eq. (20) (for small \V g \) and by Eq. (24) (for large \V g \), 
respectively. 

The ratio the transconductance g = dJ/dV g to the 
source-drain current J at high top gate voltages \V g \ > 
V b W g /W b to the same ratio at low top gate voltage \ V g | < 



FIG. 8: The source-drain current density as a function of 
the drain voltage in GNR-FETs with different gate lengths at 
different top gate voltages: solid lines - L g = 300 nm ("long" 
gate) and dashed lines - L g = 100 nm ("short" gate). 



V b W g /W b is given by 



J J high 



eV b 



low 



k B T (W b + W g ) 



> 1. 



(34) 



The source-drain current versus the back gate depen- 
dence at high top gate and drain voltages is given by 



J cx V b exp 



eV b 



k B T (W b + Wg 



1- 



1 



cosh(L 9 /2A) 



(35) 



The source-drain current increases with increasing back 
gate voltage due to the pertinent increase in the electron 
density in all regions of the channel. 



C. Short-gate effects 

As follows from the above formulas, the current-voltage 
characteristics exhibits pronounced dependence on the 
top gate length. This can be attributed to an essential 
dependence of the height of the potential barrier for the 
electrons propagating from the source on the top gate 
length (as, in particular, seen from Fig. 4) and on the 
drain voltage. As a result, the J versus Vd dependence 
(see Eqs. (29) and (31)) becomes markedly steeper when 
the top gate length decreases, particularly, when the lat- 
ter becomes comparable with A (or A) . 

In particular, as seen from Eqs. (29) and (31), in a 
GNR-FET with a long top gate, the source-drain current 
saturates when Vd becomes larger that k B T/e, whereas 
in a GNR-FET with L g comparable with A, the source- 
drain current markedly increases with increasing Vd even 
at rather large values of the latter. This is confirmed by 
the current-voltage characteristics calculated for GNR- 
FETs with a "long" (L g = 300 nm) and a "short" (L g = 
100 nm) top gates shown in Fig. 7. 
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V. DISCUSSION 
A. GNR-FET in common-drain circuit 

Above we considered the case of GNR-FETs in the 
common-source circuit (see Fig. 1), so that the voltage 
between the back gate and the source and between the 
top gate and the sourse are equal to Vb and V g , respec- 
tively. In the case of GNR-FETs in the common-drain 
circuit, the voltages Vb and V g are applied between that 
back gate and the drain and between the top gate and 
drain. In such a situation, in the above formulas one 
needs to replace Vb by V b + Vd and V g by V g + Vd- In 
this case, Eqs. (24) and (26) should be replaced by the 
following equations: 



{V g + V b W g /W b )W b 



{W b + W g ) 



V d 



1 - 



1 



cosh(L s /2A) 



V d 



2cosh(L s /2A) 



(36) 



and 



J 



2tt 3 / 2 W b 



■ exp 



k B T 



V b + V d -V b cxp[ - 



eV d 



(37) 



The main difference between the GNR-FET current- 
voltage characteristics in the common source circuit and 
in the common drain circuit is that in the latter case 
the source-drain current is a rather steep function of the 
drain voltage even at Vd S> kgT/e. Indeed, as seen from 
Eq. (31), in the GNR-FETs with L g > A, the current 
tends to saturation in the voltage range Vd 3> ksT/e. 
However, Eqs. (36) and (37) valid in the case of the com- 
mon drain circuit, give rise to an exponential voltage de- 
pendence, which at Vd ^S> ksT/e becomes as follows: 

J««p(0). (38) 



B. Role of holes 

At elevated top gate and drain voltages, the top of the 
valence band in the central section of the channel can 
become higher that the Fermi levels in the side sections. 
In this significant amount of holes can occupy the 

central section of the channel, so that a p-region can be 
formed in this section. Such an effect can occur when 
e\ip m \ > A/2 + ep or e\ip m \ + eV d > A/2 + e%. Since 
in the case of nondegenerate electron gas in the channel 
e F , e F < A/2, the condition, at which the hole gas un- 
der the top gate might be essential, can be presented as 



\<Pm\> Wm \ + Vd > A/e. As an example, assuming that 
V d < V b and setting A = 0.4 eV, W g = 0.3W 6 , and 
Vb = 2 — 3 V one obtains that the hole effects can be 
essential if \ V g \ > 1 V. It implies that the hole gas can af- 
fect the GNR-FET characteristics at elevated but realis- 
tic values of the top gate voltages. The hole charge under 
the top gate results in slowing down of the decrease of 
the source-drain current with increasing \V g \. At higher 
top gate and/or drain voltages, the tunneling between 
the n-regions in the channel not covered by the top gate 
and the p-region under this gate becomes tangible. As 
a result, the source-drain current can increase dramati- 
cally. 0,0, El BUS The tunneling current between in the 
degenerate n- and p-sections in the channel of a graphene 
FET was studied previously. @, O In the GNR-FETs 
with a marked energy band gap and a nondegenerate elec- 
tron gas in the side sections of the channel (n- regions), 
the tunneling current can become essential at fairly high 
voltages. Indeed, taking into account that the proba- 
bility of the interband tunneling between the subbands 
with n — 1 is given by 0, ll| U = exp(—n 3 hv/eE m d 2 ), 
where E m is the characteristic electric field in the n-p 
and p-n junctions arisen due to the population of the 
central section of the channel by holes. This field we 
estimate as E m = \tp m \/A = \<p m \y/3/WbW g in the n- 
p jinction and E m = (\cp m \ + Vd)y / 3/WbW g in the p-n 
junction. Using the above estimates, one can find that 
the tunneling might be essential at the values of \ip m \ 
and/or Vd exceeding the characteristic tunneling voltage 
V t = (n/4:V3)A 2 v /WbWg/ehv. Setting W b = 100 nm, 
W g = 30 nm, and A = 0.4 eV, we obtain V t ~ 6.4 V. If 
Vd <C V b and V b — 2 — 3 V, this implies that the inter- 
band tunneling might become significant at \V g \ > 9 V. 
However, the quantitative study of the tunneling effects 
in GNR-FETs with narrower energy gap, in which the 
tunneling can occur at moderate voltages, requires a sep- 
arate careful consideration which is beyond the scope of 
this paper. 



C. Main limitations of the model 

The validity of the device model used above is limited 
by the following: 

(1) Applicability of the weak nonlocality approxima- 
tion which requires that the parameter 5 — [(W^ + 
Wf)/45(W b + W g )}£ 2 < 1. This problem might arise 
only in the case of GNR-FETs with rather short top gate 
length when the terms in Eqs. (23) and (24) containing 
cosh(L s /2A) are essential. In the most interesting case, 
as follows from Eq. (22), C = ^WbW g /3. In such a case, 
S = [(Wg + W^)/15(W b + W g )]WbW g < 1. In partic- 
ular, the latter yields 5 = 1/15 ~ 0.067 if W b = W g 
and 6 ~ 0.175 if W g /Wb — 0.3. Using the Poisson equa- 
tion in the weak nonlocality approximation and retaining 
the next term in the expansion of the two-dimensional 
Poisson equation over 5, one can find that the value 
L g /2h) in Eq. (24) should be replaced by fairly close 
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value L gy /(l + 5)/2A ~ L g (l + 5/2) /2A. 

(2) Not excessively high back gate voltages (Vj, < Vp). 
This is due an assumption that the electron gas in the 
channel is nondegenerate. Setting d — 3 nm, W& = 
100 nm, ae = 4, v = 10 8 cm/s, A = 0,4 eV, and 
T = 77 — 300 K, we obtain the following estimate 
Vf — 14.0—27.6 V. This estimate, as well as the smallness 
of exp(— A/ksT), imply that inequality (10) is satisfied 
and, hence, Eqs. (14) and (15) are valid in fairly wide 
range of the back gate voltages, which definitely covers 
the range of voltages used in real devices. 

(3) Not excessively high top gate and drain voltages, so 
that the effects associated with holes under the top gate 
can be disregarded (see the previous Subsection). 

(4) The device structures with relatively week scattering 
of electrons on impurities and phonons in the active re- 
gion of the channel. The electron scattering in this region 
can give rise to a decrease in the pre-exponentional factor 
in the formula for the source-drain current in comparison 
to that in Eq. (26). However, the main voltage depen- 
dences obtained above and those which can be found us- 
ing the drift-diffusion model of the electron transport in 
the active region of the channel should be essentially the 
same. In some sense, the situation is similar to that which 
occurs in the theory of the Schottky diodes and conven- 
tional p-n- junctions when the so-called thermionic and 
diffusion models are compared [3, HH . 



One can see that the device model used above ade- 
quately describes the operation of GNR-FETs with real- 
istic parameters at reasonable applied voltages. 



VI. CONCLUSIONS 

In conclusion, we developed an analytical device model 
for GNR-FETs. The GNR-FET current-voltage charac- 
teristics, i.e., the dependencies of the source-drain cur- 
rent versus the drain voltage as well as the back gate 
and top gate voltages, were calculated for the devices 
with different geometric parameters (thicknesses of the 
gate layers and the top gate length) using this model. 
In particular, we showed that shortening of the top (con- 
trolling) gate dives rise to a substantial transformation of 
the GNR-FET current-voltage characteristics (the short- 
gate effect). The model can be used for the GNR-FET 
optimization. 
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